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Abstract 

Pion loop fluctuations of the constituent u and d quarks are shown to give 
only a minute contribution to the intrinsic charge radius of the neutron, under 
the assumption that the pion-quark coupling constant has its conventional 
value, with a cut-off scale of the order of 1.0 GeV. The contribution from 
the pion loops to the anomalous magnetic moment of the neutron represents 
a small (~ 12-14%) increase over of the static quark model value. 



The bulk of the empirical value — 0.117±0.002 fm 2 M for the mean square 
radius of the neutron is accounted for by its magnetic moment alone 0. The 
mean square radius of the neutron is defined as ||: 

<r 2 n >=-6 1imJ- 2 G E (q 2 ), (1) 

where Ge is the electric form factor. In terms of the Dirac and Pauli form 
factors G E = F 1 - (q 2 /4m 2 )F 2 , where Fi(0) = and F 2 (0) = -1.91 for the 
neutron, insertion in (1) yields the expression 

<r 2 >=-6 1imA Fl(g2) + 3^) (2 ) 
n^odq 2 2 m 2 



where m n the neutron mass. The numerical value of the latter (magnetic 
moment) term is -0.126 fm 2 , and thus is already very close to the the empir- 
ical value. The (very) small difference between this value, and the empirical 
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value, may be viewed as an intrinsic mean square charge radius of the neutron 

< r 2 > mt =< r 2 > exp -<r 2 > mag = +0.009 ± 0.002/m 2 . (3) 

This small positive value has to arise from the g-dependence of the Dirac 
form factor Fi(q 2 ). 

We here show that the lowest order pionic loop fluctuations of the con- 
stituent quarks (Fig. 1) lead to only a very small contributions to the term 
of order q 2 in Fi n (q 2 ), which yields a contribution to < r 2 >i n t, the mag- 
nitude of which with conventional parameter choices are even smaller than 
this small empirical value. 

The pion coupling to the u and d constituent quarks will be assumed to 
have the form 

£nqq = i-^^l^lpj ' (4) 

Here the pion-quark coupling constant is determined by the ttNN coupling 
constant f n NN — 1 as f^ qq ~ 3/5f n NN- As the pions should decouple from the 
constituent quarks above the chiral symmetry restoration scale A x ~ 47r/ 7r ~ 
1 GeV, a corresponding high momentum cut-off factor will be introduced at 
the pion quark vertices. 

The calculation of the pionic loop fluctuations to nucleon properties may 
be justified at the level of constituent quarks when - as in the present case 
- the smallness of the effective pion-quark quark coupling suggests a con- 
verging loop expansion. Corresponding attempts to calculate such contribu- 
tions as pion-nucleon fluctuations fail to yield realistic results because of the 
much larger pion-nucleon coupling constants and large number of intermedi- 
ate baryon resonances that have to be considered in the loops in the hadronic 
approach At the quark level the loop contributions automatically take 

into account all intermediate baryon states. 

The pionic fluctuations illustrated by the Feynman diagrams in Fig. 1 
both contribute to the q— dependent terms in the Dirac form factors of the u— 
and d— quarks. Denoting those F\ u and Fid respectively, when normalized to 
unity at q 2 = 0, the SU(Q) wave function for the nucleon yields the following 
expression for the intrinsic charge radius of the neutron: 

< r 2 n >= -4 lim ^{Fm(q 2 ) - F ld {q 2 )}. (6) 
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Figure 1: Pionic fluctuations of the constituent quarks. The fermion lines 
represent u and d quarks. Form factors are included at the vertices as ex- 
plained in the text. ^ 



Here we calculate the contributions to the Dirac form factors F lu and Fid 
from all the pion loop amplitudes that may be represented by the Feynman 
diagrams in Fig. 1, where the fermion lines represent u and d quarks, and 
the pion lines represent all the pion charge states that are allowed by charge 
conservation. Two corresponding seagull diagrams are also generated by the 
point couplings 

C^ggy = ^^((f) X f) 3 ^, (7) 

which arise by minimal substitution of the e.m. vector potential in the 
derivative coupling (4). Inclusion of the seagull diagrams is required by 
current conservation. 

In the absence of hadronic vertex form factors the seagull diagrams may be 
dropped, if in place of the derivative coupling (4) the pseudoscalar coupling: 

£nqq = ' (8) 
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is employed. The result obtained with this pseudoscalar pion-quark coupling 
is then equivalent to that obtained with the pseudovector coupling with in- 
clusion of the seagull terms for the diagrams considered here. 

This equivalence may be retained also in the presence of hadronic vertex 
form factors, provided that these are included in a way that maintains current 
conservation. Denote the vertex form factor to be inserted at each pion-quark 
vertex in the amplitudes, where the e.m. field couples to the quarks (Fig. 
lb) }'{k 2 ) (/(— m 2 ) = 1). In these amplitudes the two form factors may be 
combined with the pion propagator to a "modified" pion propagator 

<k 2 ) = -^-rr 2 f^). (9) 
mi + k 

We shall here take the form factor function f(k 2 ) to have the monopole form 
(A 2 — m 2 ) / (A 2 + k 2 ), where the value of the parameter A is ~ A x . 

To maintain current conservation, and the equivalence between the pseu- 
dovector (including seagull terms) and pseudoscalar coupling models, the 
vertex factors have to be inserted into the amplitude that corresponds to the 
diagram Fig. la, where the e.m. field couples to the pion current, so that the 
product of the two pion propagators is modified to |§, 0] : 

1 1 : v(k 2 )-v(k 2 ) 

ml + klml + kl k\-k\ ' 

Here k\ and k 2 represent the 4-momenta of the pion before and after the 
electromagnetic coupling. Form factors inserted by this method allows the 
calculation to proceed on the basis of the pseudoscalar coupling (8) without 
further reference to seagulls. 

The e.m. vertex of the internal constituent quarks, which should be 
assumed to have a spatially extended structure that may be described by a 
form factor F q (q 2 ), is described by the current matrix element 

< PU(«, d)\p >= ze(| -l)u(p'){^ + (F q (q 2 ) - 1)[ 7 „ - 1 -^}}u(p). (11) 

This vertex, in which the first term corresponds to pointlike quarks, main- 
tains the requirement of current conservation, as the form factor modification 
appears only in a purely transverse term [|, ||. 



4 



The pion should similarly be described as spatially extended, in view of its 
empirically large mean square radius. This is well described by the p— meson 
pole in the time-like region. This is taken into account through inclusion of 
the pion form factor F n (q 2 ) in the e.m. coupling of the pions, which is taken 
to have the form 

< k'\j^°)\k >= ie(±l,0){K, + (F^q 2 ) - l)[K, - ^f%, (12) 

where = + In this vertex the form factor also appears only in a 
transverse term. Note as the electromagnetic vertices are considered at tree 
level, they may contain all information of the quark structure of the pion 
and the self-dressings of the pion and the constituent quarks without double 
counting. 

Consider the pionic fluctuations of u and d quarks illustrated in Figs, la 
and lb. The only processes, in which the the e.m. field couples to the pions, 
are the fluctuations u — > n + d — > u and d — > tt~u — > d. The flavor factors of 
these two processes are the same, and therefore their contributions to the u 
and d quark current matrix elements have the same magnitude, but opposite 
sign. 

For u quarks the fluctuations with e.m. coupling to the internal quark 
lines are u — > 7i + d — > u and u — > ii u — > u. In the first the coupling 
is proportional to the d quark charge — e/3, while the flavor factor from 
the quark charge changing vertices is (\/2) 2 . In the second the coupling is 
proportional to the u quark charge 2e/3, while the flavor factor from the 
quark charge conserving vertices is 1. Multiplication of the charge and flavor 
factors shows that these two contributions cancel exactly in the case of the 
u quark. In the case of the d quark the product of charge and flavor factors 
for the corresponding fluctuations d — ■> n~u — > d and d — > n°d — > d are on 
the other hand 4/3 and —1/3 respectively, which add up to 1. Consequently 
fluctuations with e.m. coupling to the intermediate quark only contribute to 
the current matrix element of the d quark. 

The evaluation of the net loop contributions to the form factor combina- 
tion Fiu — Fu required for the neutron mean square radius (6) is particularly 
convenient, as the combined current matrix elements satisfy the continuity 
equation, without e.m. coupling to the external quark legs. In the absence 
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of form factors one obtains the result: 



^ r n ^loops— o I / dxX { . . 2 



7r 2 I Jo H{mir) H{m 1T ) 2 

Here the function H(m) is defined as 

H(m) = m 2 q (l - x) 2 + m 2 x. (14) 

The first integral on the rhs of (13) arises from fluctuations with e.m. coupling 
to the pions (Fig. la) and the second from terms with e.m. coupling to the 
quark lines (Fig. lb). The latter is finite in the chiral limit, while the former 
is divergent in that limit. That divergence also appears in other chiral pion 
field theoretical models, e.g. in the Skyrme model ||. 

The monopole hadronic form factors at the pion-quark vertices introduced 
so as to maintain current conservation as described above are included by 
means of the following substitutions in the expression (13): 

1 (16.) 



H(m w ) H(jn n ) H(A) H(A) 2 

1 1 1 A 2 — m 2 

2x TT ™* . (156) 



ff(m,) 2 H(m^) 2 H(A) 2 if (A) 3 

The pionic loop contribution to the mean square neutron radius (13) is 
very small and sensitive to the constituent mass value. With a quark mass 
value of m q = 300 MeV the numerical values are -0.00082 fm 2 and 0.0018 
fm 2 for A = 800 MeV and 1.0 GeV respectively. With m = 340 MeV the 
corresponding numerical values are -0.0048 fm 2 with A = 800 MeV and - 
0.0023 fm 2 with A = 1 GeV. The value is positive for A > 1.2 GeV, the 
asymptotic value for A — >• oo being 0.013fm 2 with the latter constituent mass 
value. This latter value is close to the empirical range. The smallness of the 
pion loop contribution to the mean square neutron radius is a consequence 
of strong cancellations between the two terms in (13). 

The additional contribution to the mean square radius of the neutron 
from the finite radii of the pion and the quarks within the loops are 

9 <?,9laas r 1 , fir, H(A) A 2 -TO 2 
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-^-'^H^rrm-'^w^ (16o) 



a/ sLv f 1 , vfl n #(A) A 2 -m 2 

A 2 - m 2 



^-^fe-tf^W*)' (166) 

respectively. 

The mean square pion radius obtained from the empirical parameters 
in ref.[]lU) is r 2 = 0.38 fm 2 . Realistic dynamical models for the baryon 



spectrum |TT suggest that the (flavor averaged) mean square matter radius 
of the constituent u and d quarks should be about 0.13 fm 2 if the empirical 
mean square radius of the proton is to be reached JT^]. With these values 



the additional contribution to the mean square radius of the neutron from 
the pion and (flavor averaged) constituent quark radius is very small. With 
m q = 300 MeV the value is 0.00086 fm 2 for A = 800 MeV and 0.0020 
fm 2 for A=1.0 GeV. With m q = 340 MeV the corresponding values are 
somewhat smaller: 0.00013 fm 2 and 0.0014 fm 2 for A = 800 MeV and A = 
1 GeV respectively. The smallness of these values is again a consequence of 
cancellations between the two contributions (16a) and (16b). 

Combination of the two contributions (13) and (16) to the neutron radius 
then gives the following net contributions to the mean square radius of the 
neutron: with the quark mass 300 MeV and A = 800 MeV the net contribu- 
tion is only 0.00004 fm 2 , whereas with A = 1 GeV it amounts to 0.004 fm 2 . 
In the case of m q = 340 MeV the net contribution is -0.005 fm 2 for A = 800 
MeV and -O.OOlOfm 2 for A = 1 GeV. With the smaller quark mass value the 
net contribution reaches 0.00 fm 2 and thus the empirical range 0.009±0.002 
fm 2 if the cut-off parameter A is increased to 1.2 GeV. With the larger quark 
mass value the empirical range is reached only by A = 1.5 GeV. The conclu- 
sion is in any case that the pionic loop fluctuations of the constituent quarks 
imply only a very small contribution to the intrinsic mean square radius of 
the neutron. The pionic loop contribution does therefore not perturb the 
satisfactory description of the (bulk of the) negative mean square radius of 
the neutron, which is implied by the empirical value of the neutron magnetic 
moment. 

The pionic fluctuations of the constituent quarks also give a small contri- 
bution to the magnetic moment of the neutron. In units of nuclear magnetons 
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this contribution is 



g 2 f 1 1 1 A 2 — m 2 

F 2 (0) loops = -^rn p m q J o Ml-^+SxKj^-j^-x^^}. 

(17) 

Here m p is the proton mass. The numerical value of this contribution is only 

0.23 n.m. for A = 800 MeV and 0.26 n.m. for A = 1 GeV with 

m q = 340 MeV. The asymptotic value for A — > oo is 0.37 n.m. These values 
represent enhancements of about ~ 12-14% of the static quark model value 
—2/3(m p /m q ). The latter agrees with the empirical value -1.91 n.m. with 
m q = 345 MeV, but only if the reduction caused by the considerable "rel- 
ativistic correction" that appears if the quarks are described by the Dirac 
current operator is neglected |15j . Hence there is a room for this loop correc- 
tion along with exchange current contributions that are associated with the 
interaction between the confined quarks |T2, |T^, [TJ[ . 

The pion loop contribution to the anomalous magnetic moment of the 
proton -Ff (0) may also be calculated using the expression (17), provided that 
the bracket (2+3x) in the integrand is replaced by the expression -(l+9x)/2. 
Using the same parameter values as above for the neutron, we find the loop 
contribution to F$ (0) to fall in the range 0.17 - 0.19 n.m. The static quark 
model value for the magnetic moment of the proton is m p /m q , and for the 
ratio of the neutron to proton magnetic moments = 2/3. This ratio, 

which is slightly below the empirical value 0.68, would increase by ~ 5% by 
inclusion of the loop contributions considered here. 

The pionic loop contribution to the neutron magnetic moment was re- 
ported to be much larger than the value above in refs. |TB|, [T7]], without 
supporting formalism. The suggestion in refs. [[16], |17| that the contribution 
of the pionic (and kaonic and rj— meson) loop fluctuations of the constituent 
quarks to their anomalous magnetic moments explain all of the negative em- 
pirical mean square radius of the neutron is however not tenable as it neglects 
the other (main) quark contributions that make up the (empirical) anoma- 
lous magnetic moment of the neutron. In addition the contribution to the 
neutron mean square radius from the Dirac form factors of the constituent 
quarks to the first term of eq. (2) was neglected in refs. |L6|, fL7|| . 

The present finding is that the lowest order pionic loop fluctuations are 
insignificantly small for for conventional values of the coupling and cut-of 
scale parameters. These loop corrections are the leading ones in a 1/m 
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expansion, being of order (1/m) 2 . Small additional contributions that are of 
order (1/m) 3 arise from exchange currents [13, 14 . 



The restriction to SU (2) flavor symmetry above is not essential, and the 
Ks and rjs fluctuations of the u and d quarks may calculated by methods 
similar to those used above. The much larger masses of the K and r\ mesons 
make the numerical contributions of those fluctuations smaller than those of 
the pionic fluctuations however. Vector meson loop fluctuations are expected 
to be small, because of the much larger meson masses. Loop fluctuations in 
which intermediate vector mesons undergo radiative decay to pseudoscalar 
mesons give no contribution at all to the intrinsic charge radius of the neu- 
tron. This is direct consequence of the transversality of the transition current 
matrix elements, which have the generic form 

< ir(k')\ -h\V h M >= i 9 -^e, Xua k x k' u 6 ab . (18) 

Here V£ denotes either the p or the uo meson field (in the latter case the 
isospin index is left out). 
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